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(J , We asymptotically compute the intersection angles between A^- dimensional stable 

and unstable manifolds in 2A^-dimensional symplectic mappings. There exist particular 

/\ • 1-dimensional stable and unstable sub-manifolds which experience exponentially small 

c^ I splitting of separatrix in our models. We show that the angle between the sub- manifolds 

is exponentially small with respect to the perturbation parameter e, and the other 

angles are O(e^). 
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The existence of a transversal intersection between stable and unstable manifolds leads to 
non-integrability of a Hamiltonian system. Besides, in 2-dimensional symplectic mappings 
(area-preserving mappings), which correspond to Poincare mappings of Hamiltonian systems 
with 2 degrees of freedom, the area of the lobes enclosed by stable and unstable manifolds 
stands for flux from inside (resp. outside) of separatrix to outside (resp. inside). Because 
the intersection angle at the principal intersecting point is related with the area of the lobes, 
the intersection angle is also related with flux. Therefore the intersection angle is essentially 
important to understand flux near separatrix and, probably, global transport of phase space. 

For some 2-dimensional symplectic maps, e.g. the standard map, Henon map, and the 
double-well map, one can analytically construct asymptotic expansions of stable and unstable 
manifolds. a^'cP'U^'H' Hence one can also compute the asymptotic form of the intersection 
angles between stable and unstable manifolds in such maps. We have been extending the 
method to 4- or more- dimensional symplectic maps with two or more hyperbolic modes, to 
which Melnikov's methodH-* cannot be applied. 

Such an approach was first presented by Gelfreich and SharomovQ^ and we generalized 
their model in Ref. |^). In addition, we constructed the asymptotic expansions of partic- 
ular 1-dimensional stable and unstable sub-manifolds which experience exponentially small 
splitting of separatrix, in a 4-dimensional double-well symplectic mapping. Furthermore, the 
neighborhood of the sub-manifolds was asymptotically studied and exponentially small oscil- 
lating terms were successfully obtained. eP In this paper we treat 2A^-dimensional double-well 
symplectic mappings and standard mappings with nearest-neighbor two-body couplings and 
asymptotically compute the intersection angles between iV-dimensional stable and unstable 
manifolds. 

The canonical coordinates we use are q= (gi, ■ ■ ■ , Qn) ^ R^, P = {Pir ' ' ^Pn) ^ R-^ and 
7 = (^)P)- We start with 2iV-dimensional symplectic mappings 7 1— >• 7', 

^^ - ^^ % % ' (1) 

q'j = Qj + ep'j, 

where j = 1,2,- ■ ■ , N^ and |e| ^ 1 is a perturbation parameter. We take account of 
two different mappings: One is the double- well mapping (DW), the other is the standard 

mapping (SM). The potential function Vo(gj) is given as -q'^ g^ (for DW), cosg (for SM), 

respectively. For both of the potential functions, the origin (0, ■ ■ ■ , 0) is a 2A^-dimensional 
hyperbolic fixed point (a direct product of A^ 2-dimensional hyperbolic fixed points), and 
there exist A^-dimensional stable manifolds and A^-dimensional unstable manifolds near the 



^ We shall use the subscript j as meaning j = 1,2,- ■ ■ ,N without further notice. 



origin. The coupling potential Vi{q) is an even polynomial of the distance between nearest- 

^ a- ^ b- 

neighbors of order 4, i.e., Vi{q) = ^ i^ilj ~ Qj+iY + X! ^(^i ~ 9'i+i)''' where a/s and 6j's 

are real constants and g^v+i = Qi- In this paper we restrict ai = a2 = ■ ■ ■ = ajy = a and 
hi = 1)2 = ■ ■ ■ = h^ = h to simplify the following computation. 

Setting qj{t) = Qj and qj(t + e) = q'j, one can rewrite the map (1) into the second order 
difference equation 

A',q,{t) = U{q,) + e'l)M (2) 



T-^-'^-^-l" 



where Z\2^(t) = {q{t + e)-2q{t) + q{t-e)}/e^ /o(g,) = £^ and /,, i(g-) = ^. By 

regarding the independent variable t as continuous one, the difference equation (2) expresses 

time evolution of not trajectories but invariant manifolds, because the set of trajectories 

°° {±ty d}q- 
between t and t+e can be took into consideration in (2). Expanding gj(t±e) = ^ — — y-, 

one obtains 

^^/»fe)+eV..®-2|:L_^. (3, 

which we call the outer equation. 

The stable and unstable solutions q^ (t; e) are asymptotic to the origin as t ^ ±oo, 
respectively. Hence the boundary condition with which we solve the difference equation (2) 
is given as qf{t; e) —^ +0 (for DW), =fO (for SM) as t tends to ±oo. 

We construct the stable and unstable solution to the outer equation (3) by expanding 
qf{t] e) into the power series of e^, i.e., qf{t^ e) = qf^it) + e^qfi{t) + ■ ■ ■. In the same way, 
the stable and unstable manifolds 7^ is expanded as 

7^ = 7o'= + e'7i^ + ---. (4) 

By substituting the expansion into the outer equation (3), one gets the unperturbed 
equation as 

# = /»(«|i>) (5) 

and the perturbed equation (the linearized equation) as 



^ = mMi + hAt) (6) 



for / = 1,2, ■■■, where the inhomogeneous term fj^i is a polynomial of qj^Q{t)^s and their 

^, 1 d^qjQ 
derivatives, e.g., fj^i = fj,i{%) 74 • ^"^ ^'^^ same way the boundary condition is 



also expanded as qfi -^ +0 (for DW), =fO (for SM) as t ^ ±00, / = 0, 1, ■ ■ ■. Note that 
the ordinary differential equations for g^y, / = 0, 1, ■ ■ ■ form a Hamiltonian system with N 
degrees of freedom in order-by-order. 

The unperturbed solution qf^it) is obtained as q^Q = s{t + tj) where s{t) = sech(t) (for 
DW), 4tan^-'^(e*) (for SM), and tj's are constants of integration. Let us put sit + tj) = Sj(t) 
to simplify notations. The unperturbed solution Sj(t) does not depend on the sign ±, i.e. 
the stable and unstable solutions behave completely same for —00 < t < 00 and intersect 
tangently. That is a manifestation of integrability of the unperturbed equation (5). 

By taking A^ of tj's and t as independent parameters, A^-dimensional stable and unstable 
manifolds can be constructed. Because the map (1) is autonomous, it is translationally 
invariant with respect to t. Hence one of t^'s is not essential and can be eliminated. This 
point will be argued later. 

To construct the solutions to the linearized equation (6), we first consider the homo- 

(Pqfi , 

geneous equations, — ry- = foicijflJQj,^ which are decoupled into A^ independent equations. 

The fundamental system of solutions to the homogeneous equation is Sj{t) = s{t + tj) and 

f dt 
Qjit) = g{t + tj) where g{t) = s i ^ . The stable and unstable solutions to the equation (6) 

can be written as 

rt rt 



&^ti) = 9j{t) / s,(t)4i(t)rft - s,{t) / gj{t)f,^i{t)dt, (7) 

J±oo Jo 

respectively. 

When tj = 0, the r.h.s. of the expression (7) can be easily integrated. The result 
does not depend on j and the sign ±, hence putting gi(t) = qf-^^(t,0), one obtains qi{t) = 

-sech^(t) - — sech(t) + — sinh(t)sech^(t) (for DW), -sinh(t)sech2(t) - — sech(t) (for SM). 
The solution qf{t, 0) = s{t) + e'^qi{t) + 0{e^) stands for the 1-dimensional stable and unstable 



^3 

sub-manifolds, which we call separatrix. Note that the sub-manifolds experience exponen- 
tially small splitting and split into two different manifolds in the order of exp .EP 

To simplify the further computation we transform t^-'s into new parameters a^-'s as 
(ti, • • • , t^) = («!, ■ ■ ■ , aN)M, where 

f -^ 1 \ 

-1 ■•. 
■•• 1 

Because detM = (—1)^ 7^ 0, the transformation is non-degenerate. We put aN = without 



M 



loss of generality. The A^-dimensional stable and unstable manifolds are spanned with these 
A^ — 1 parameters, ai, ■ ■ ■ , aN-i and t. 

The tangent spaces of the sub-manifolds are necessary to estimate the intersection an- 
gles on the sub-manifolds. We write the stable and unstable manifolds as 7^(t,aj) = 
(g^(t, (Jj), ^^^(t, (Jj)), respectively. Let us suppose that the stable and unstable manifolds 
intersect at a point zq = 7^ (to, 0) where to G R. We write tangent vectors of the manifolds 

at the point as Xf, where X^ = -^{to, 0), fc = 1, ■ ■ ■ , iV - 1 and X^ = ^— (to, 0). We 

— * fi IT 

write X^ = (a^(to),'i^(to)), where v^{t) = —r~, and x^{t) can be expanded with respect 

to e as f^(t) = x^fc(t) + e^x^f,(t) + O(e^). Owing to 70 = (si(t), ■ ■ ■ , ■SAr(t), si(t), ■ ■ ■ , s„(t)) 

(see the expression (4) and the discussion below it), one obtains 

d 

Xo,k = ^ (si(t),---,SAr(t)) 

dak 5=0 

k fc+1 

= (0, • • ■ , 0, -s, s , 0, ■ ■ ■ , 0), for A; = 1, ■ ■ ■ , A^ - 1 
Xo,N = {s,---,s)- (9) 

Note that X^ ± X^, X^ ± X^+„, A; = 1, ■ ■ ■ , A^ - 1, m > 2, and X^ / X^+i, k^ N -1 
up to 0(e2). 

The set of tangent vectors {X • } is a basis of the tangent space of the stable and unstable 
manifolds at zq, respectively. We denote the inner product between X+ and X~ as (X+|X~) 
where 1 < m,n < N. The intersection angles between the stable and unstable manifolds at 
Zq are defined as 

emAto) = COS-l^S^. (10) 

Iml In I 

Especially, 9]\j^n stands for exponentially small splitting of separatrix,a''}3'Q> i.e., there exist 
C > such that 

\ON,Nito)\<^exp(-^y (11) 

where a = 5 (for DW), 3 (for SM), at the principal intersecting point (to = 0). 

To compute the intersection angles 9m,n, one has to estimate X^ including the terms of, 
at least, O(e^). Hence we shall analyze the neighborhood of the separatrix as a perturbation 
problem with respect to a/s. We expand qfA^,a) (see the expression (7)) with respect to 



aj s as 



N-l 

E 

fc=0 



g± (t,a) = gi(t) + ^ a^qf^t) + Oi\a\''), (12) 



where qfk{t) = ^' 



Oik 



. From the expression (7), straightforward computation gives 



a=0 



q%{t, a) = t,,kFt{t) + (t,+i,fc + t,.,,k)F^{t), (13) 

da} 
s{t) and g{t) as 



where tj^ = -^^- The functions F^ [t] and Fg (t) is exphcitly written by using the functions 
oak 



F^{t) = qi{t) ± aPg{t) - a5s{t) sinh2(t), 

F^t) = T\aPgit), (14) 

/oo 2 1 

s{tfdt = - (DW), 8 (SM), and (5 = - (DW), 1 (SM). The detail of this com- 
-oo O I — [ O 

putation will be reported elsewhere. liHl' Note that the expression (14) contains the parameter 
a, which determines the strength of coupling. 

We are now ready to give the explicit representation of Xf" including the terms of O(e^). 

Because 71 = (qf^i, ■ ■ ■ , qN,i^Pf,ii ' ' ' iPn,i) where pf^^ = J' , by using the expression (13), 
one can explicitly write xf,^ as 



^l,k — (Q'i,15 ■ ■ ■ 5 Qn,!) 
k-l 

= {0,---,0,-Ft,F^-F,^,F,^-F^,F^,0,---,0), (15) 



fc-l k k+1 fe+2 



for k = 1,- ■ ■ , N — IM By substituting the expressions (9) and (15) into the equation (10) 
and using the expression (14), one obtains the intersection angles between the stable and 
unstable manifolds at the point zq as 

OkAto) = .,, ,ff ■' , ,, [9 + (^(to)^ + 9{tor){s{tor + s{tor)f^ + 0{e'), 

S[to) + S[tQ) 

ek,k+i{to) = Y + 0(e2), for A; ^ iV - 1 
Ok,k+m{to) = ^ + O(e'), for m > 2 
OkAto) = l + 0{e'), (16) 

for fc = l,---,iV-l. 

^ The index of the elements is taken periodically and overlapped elements (if exist) are summed. For 
example, when N — 2, x^^ = (— -F^ + 2F2 , Ff^ — 2F2 ) by definition and the expression (13). The first 
element is regarded as the sum of the elements of index k and fc + 2(fc = fc + 2 mod 2) , and in the same way, 
the second element is the sum of index fc — 1 and fc + 1. 



The angles between the neighbors are not orthogonal with each other. That is based on 
the fact that Xk is not orthogonal to Xk+i, and is not essential. Although one can let them 
orthogonal by changing the matrix (8) into others, the notations will become messy. 

The expression (11) and (16) imply that A^ of the angles between the stable and unstable 
manifolds, i.e. Ojj, j = 1,---,N are independent and sufficient to write down the inter- 
section of these A^-dimensional manifolds, and the results tell one that the only one angle 
is exponentially small and the others are power of e. In the outer approximation (an ap- 
proximation of the difference equation (2) by the the differential equation with the singular 
perturbations (3)), the outer equation (3) forms a Hamiltonian system with A^ degrees of 
freedom in order-by-order and possesses an "energy function." Therefore trajectories of the 
outer equation are restricted in {2N — l)-dimensional energy surface, and in the energy sur- 
face the intersection of the A^-dimensional stable and unstable manifolds has 1-dimension. 
In reality, because we consider not continuous Hamiltonian systems but symplectic map- 
pings, then such energy functions do not exist and the separatrix splits into two different 
manifolds, i.e., the 1-dimensional stable and unstable manifold, in the exponentially small 
order. EP We think that this feature is quite generic over the symplectic mappings. The naive 
approximation misses the exponentially small terms and leads to an incorrect result that the 
intersection angle between sub-manifolds is equal to zero. The approximation is however 
sufficient to estimate the angles of the directions out of the separatrix. 

Intersection angles are related with flux near stable and unstable manifolds. We expect 
that escape rate of trajectories near the manifolds gets large as A^ increases because the 
number of the angles in the power of e dominates overwhelmingly. We would like to conflrm 
this conjecture with numerical experiments in the future work. This point is closely related 
with well-known Nekhoroshev's method and a very important problem for understanding the 
global structure of phase space of nearly integrable Hamiltonian systems with many degrees 
of freedom. 
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